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Abstract: Exclusive semileptonic B decays to orbitally and radially excited charmed mesons are 
investigated in the first order of the heavy quark expansion. The emerging leading and subleading 
Isgur-Wise functions are calculated in the framework of the relativistic quark model. It is found that 
both relativistic and the l/mg corrections play an important role and substantially modify results. 
An interesting interplay between different corrections is observed. 
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1. Introduction 

The investigation of semileptonic decays of B 
mesons to excited charmed mesons represents a 
problem interesting both from the experimental 
and theoretical point of view. The current exper- 
imental data on semileptonic B decays to ground 
state D mesons indicate that a substantial part 
(« 40%) of the inclusive semileptonic B decays 
should go to excited D meson states. First exper- 
imental data on some exclusive B decay channels 
to excited charmed mesons are becoming avail- 
able now H H and more data are expected 
in near future. Thus the comprehensive theoret- 
ical study of these decays is necessary. The pres- 
ence of the heavy quark in the initial and final 
meson states in these decays considerably simpli- 
fies their theoretical description. A good start- 
ing point for this analysis is the infinitely heavy 
quark limit, tuq — ► oo 0. In this limit the heavy 
quark symmetry arises, which strongly reduces 
the number of independent weak form factors jf| . 
The heavy quark mass and spin then decouple 
and all meson properties are determined by light- 
quark degrees of freedom alone. This leads to a 
considerable reduction of the number of indepen- 
dent form factors which are necessary for the de- 
scription of heavy-to-heavy semileptonic decays. 
For example, in this limit only one form factor 
is necessary for the semileptonic B decay to S- 
wave D mesons (both for the ground state and its 



radial excitations), while the decays to P states 
require two form factors ||. It is important to 
note that in the infinitely heavy quark limit ma- 
trix elements between a B meson and an excited 
D meson should vanish at the point of zero re- 
coil of the final excited charmed meson in the rest 
frame of the B meson. In the case of B decays 
to radially excited charmed mesons this follows 
from the orthogonality of radial parts of wave 
functions, while for the decays to orbital excita- 
tions this is the consequence of orthogonality of 
their angular parts. However, some of the l/mg 
corrections to these decay matrix elements can 
give nonzero contributions at zero recoil. As a 
result the role of these corrections could be con- 
siderably enhanced, since the kinematical range 
for B decays to excited states is a rather small 
region around zero recoil. 

Our relativistic quark model is based on the 
quasipotential approach in quantum field the- 
ory with a specific choice of the quark-antiquark 
interaction potential. It provides a consistent 
scheme for the calculation of all relativistic cor- 
rections at a given v 2 /c 2 order and allows for the 
heavy quark 1/rriQ expansion. In preceding pa- 
pers we applied this model to the calculation of 
the mass spectra of orbitally and radially excited 
states of heavy- light mesons M, as well as to a 
description of weak decays of B mesons to ground 
state heavy and light mesons |t], |j. The heavy 
quark expansion for the ground state heavy-to- 
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heavy semileptonic transitions ]9| was found to 
be in agreement with model-independent predic- 
tions of the heavy quark effective theory (HQET). 

2. Decay matrix elements and the 
heavy quark expansion 

In this section we present the heavy quark expan- 
sion for weak decay matrix elements between a B 
meson and radially excited charmed meson states 
up to the first order in 1/toq using the HQET. 
The corresponding formulars for B decays to or- 
bital excitations can be found in Ref. |]Io| . 

The matrix elements of the vector and axial 
vector currents between B and radially excited 
D' or D*' mesons can be parameterized by six 
hadronic form factors: 



(D'(v')\c^b\B(v)} 



h+(v + v'y + h-{v- v'Y, 



(D'(v')\c^ l5 b\B(v)) = 0, 



y'lfiD-'nifl 
{D*'{v',e)\c^ l5 b\B{v)) 



h Al (w + l)e*» 

{h A2 v» + h A3 v'»)(e* -v), 



(2.1) 



where v (v') is the four-velocity of the B (i)W') 
meson, e M is a polarization vector of the final vec- 
tor charmed meson, and the form factors hi are 
dimensionless functions of the product of veloci- 
ties w = v ■ v' . 

Now we expand the form factors hi in powers 
of 1/rriQ up to first order and relate the coeffi- 
cients in this expansion to universal Isgur-Wise 
functions. This is achieved by evaluating the ma- 
trix elements of the effective current operators 
arising from the HQET expansion of the weak 
currents. For simplicity we limit our analysis to 
the leading order in a s and use the trace formal- 
ism fll]. Following Ref. go|, we introduce the 
matrix 



H,, = 



i+ i> 



P„*" 7m - P„ 7 5 



(2.2) 



composed from the fields P v and P* M that de- 
stroy mesons in the j p = \ doublet 1 with 
four- velocity v. At leading order of the heavy 



iHere j is the total light quark angular momentum, 
and the superscript P denotes the meson parity. 



quark expansion (jtiq — > oo) the matrix elements 
of the weak current between the ground and ra- 
dially excited states destroyed by the fields in H v 
and H' v , respectively, are given by 



cTb 



h$ThV> =^\w)Tr{H^TH v }, (2.3) 



where is the heavy quark field in the effec- 
tive theory. The leading order Isgur-Wise func- 
tion £( n >(w) vanishes at the zero recoil (w = 1) of 
the final meson for any T, because of the heavy 
quark symmetry and the orthogonality of the ra- 
dially excited state wave function with respect to 
the ground state one. 

At first order of the l/mq expansion there 
are contributions from the corrections to the 
HQET Lagrangian 



SC 



1 



(Q) 



1 



2mc. 
q(Q) . 

Q(Q) 
v mag,ii 



q(Q) ,q(Q) 



" 2m Q 



,(2.4) 



and from the tree- level matching of the weak cur- 
rent operator onto effective theory which con- 
tains a covariant derivative D x = d x — ig s t a A x 



cTb 



r 



2m, 



■pT- 



2rrifc 



r^Uw. (2.5) 



The matrix elements of the latter operators can 
be parameterized as 



(2.6) 



The most general form for £ x is [O 



AQ) 



e 



(Q) i 



'(v + v^x+e^iv-v')^-^'^- 

(2.7) 

The equation of motion for the heavy quark, 
i(v ■ D)h^ — 0, yields the relations between the 
form factors £^ 

ei c) (i + H+d c) (w-i)+4 c) = ° 

^\l + w)-e\w-l) + 4 b) = 0. (2.8) 

The additional relations can be obtained from 
the momentum conservation and the definition 
of the heavy quark fields hi^ , which lead to the 
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equation id v {h[ c )v h^ ] ) = (Au !/ -A(")^)/i^ ) r^ b) , 
implying that 



(2.9) 

Here A(V\>)) = M(M^) - m Q is the differ- 
ence between the heavy ground state (radially 
excited) meson and heavy quark masses in the 
limit rriQ — > oo. This equation results in the fol- 
lowing relations 



(b) 



.Ac) 



Ac) 



A b ) _ A c ) _ ^ 



Ac) 

S3 



d b) = A£ (n) , 

(&) = __^(n)^(n) ) 

= 0. (2.10) 



Ab) 
S3 



The relations (2.8) and (2.10) can be used to ex- 
press the functions £_ 1 in terms of £3 (= £3 ^ = 



£3 ) and the leading order function £ 



(n). 



.(c) 



^(b) = 

<r(c) _ 
s+ — 



A(«) 1 A( 

~Y~ + 2 w- 1 

A 1 A(") - A 
2 ~~ 2 



A 



A(«) 



w — 1 
1 AW 

f 2 — 



A 



(6) 



1A<") 
2 w- 



w; + 1 
-A 



:(») 



r6, 



1 



(n) 



1 



6- (2.11) 



The matrix elements of the 1 / tuq corrections 
resulting from insertions of higher-dimension op- 



erators of the HQET Lagrangian (2.4) have the 
structure Q 

dxT{4% (x) [h&rW] (0)} 
= 2 X l c) Tr{^r^} 



+2Tr | X giW a/5 

d xT {c[%) [TtfrW] (0)} 

= 2 X f ) Tr{H v ,TH v } 

+2Tr \ X %B vl T^-lia^H t 



,1+ / 



(2.12) 



The corrections coming from the kinetic energy 
term Okin do not violate spin symmetry and, 
hence, the corresponding functions X i effec- 
tively correct the leading order function ^ n > . The 
chromomagnetic operator O mag , on the other 
hand, explicitly violates spin symmetry. The most 
general decomposition of the tensor form factor 



X a /3 



(c) 

X2 V alf3 



(c) ■ 
X 3 



(*>) (b) 1 (b) . in ^1\ 

Xjp = X2 v aJ0 ~ X3 ™<x$- (2.13) 



The functions Xi contribute to the decay 



form factors (2.1) only in the linear combination 
Xb = 2xf ) -4(w-l)x2 fc) + 12x^ ) . Thus five inde- 
pendent functions £3, Xfc and x^(= Xi^)i as we ^ 
as two mass parameters A and A' n ) are neces- 
sary to describe first order 1/toq corrections to 
matrix elements of B meson decays to radially 
excited D meson states. The resulting structure 
of the decay form factors is 

h+ = + e c [2xi - 4(u> - 1) X 2 + 12X3] + e bXb , 



26 



A («) + A(n) Z A ) £(") 
w — 1 



-E6 

-4X3 



U> — 1 

2 A i I A< 



AW - A 



1/7 — 1 



i(n) 



w — 1 



A 



A(«) - A 



w — 1 



26 



2xi - 4X3 



+ 



+£& i Xb + 



w + 1 \ iu — 1 

u> — 1 



10 + 1 
/lA 2 = £c<{ 4x2 - 

+6 



a-^AV«> 

w - 1 ' 



26 



U) + 1 



A(n) + A(W) ~A \ £ (n) 
w - 1 ' " 



10 + 1 



+£ b 



Xb 



A 



2xi - 4x2 - 4x3 

" AW - A s 

w — 1 
A(«) - A 



£(«) _ 



U) + 1 



6 



w — 1 



26 



(2-14) 



where sq = 1/(2?tiq). 

The similar analysis |l(| for B decays to or- 
bitally excited states indicate that it is necessary 
to introduce two Isgur-Wise functions in leading 
order of the heavy quark expansion: one function 
t(w) for decays to Di^D^ mesons with j = 3/2 
and the other one C( w ) f° r decays to Dq,DI 
mesons with j = 1/2. At subleading order six 
additional functions (r\ t 2, rjk e , ^1,2,3) arise for the 
former decays and four functions (£1, Xfce, Xi,2) 
for the latter ones. 
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3. Relativistic quark model 

We use the relativistic quark model based on 
the quasipotentil approach for the calculation of 
corresponding Isgur-Wise functions. Our model 
has been described in detail at this conference 
H|, so we directly go to the calculation of de- 
cay matrix elements of the weak current between 
meson states. In the quasipotential approach, 
the matrix element of the weak current J = 
67^(1 — 7 5 )& between a B meson and an excited 
D** meson takes the form 

(D**\JW(0)\B) 

^^ D ,.(p)r M (p,q)vMq), (3.1) 

where T AI (p, q) is the two-particle vertex function 
and ^b.d** are the meson wave functions pro- 
jected onto the positive energy states of quarks 
and boosted to the moving reference frame. The 
contributions to T come from Figs. 1 and 2. 2 In 
the heavy quark limit m(, jC — > 00 only con- 
tributes, while r' 2 ) contributes at 1/mg order. 
They look like 

rW(p,q) = u c ( Pc )^(l- 1 5 )u b (q b )(27:) 3 d(p q -c iq ), 

(3.2) 

and 

r^ 2) (p,q) = u c {Pc)u q {p q ){~f Q v(l - 7q) 



eb(k) + e b (p c 
e c {k') + e c {qb) 



-7qV(p 9 - q q ) + V(p q - q q ) 

7q7Qm(! - jQ)} u b(Qb)uq(q q ), 

(3.3) 



where the superscripts "(1)" and "(2)" correspond 
to Figs. 1 and 2, Q = c or b, k = p c — A; k' = 
q h + A; A = p D „ - p B ; e(p) = (to 2 + p 2 ) 1/2 ; 



A(-)(p) = 
Here [H 



i(p) - (to7° + 7°(7P)) 

Hp) 



Peg = e c , q (p)j^- ±Y^n ii) (p D *>)p i 



i=l 



2 The contribution r' 2 ) is the consequence of the pro- 
jection onto the positive-energy states. Note that the 
form of the relativistic corrections resulting from the ver- 
tex function r( 2 ) is explicitly dependent on the Lorentz 
structure of the gq-interaction. 



W 



B 



D* 



Figure 1: Lowest order vertex function r' 1 ' cont- 
ributing to the current matrix element (3.1). 



•J 

B 



W 



D* 



B 



j W 

D* 



Figure 2: Vertex function F' 2 ' taking the quark 
interaction into account. Dashed lines correspond 
to the effective potential. Bold lines denote the ne- 
gative-energy part of the quark propagator. 



qb.q = e b , q {q)^f-±y^n {l) {p B )q\ 
i—i 

and are three four-vectors given by 



»«"(p) = 



M' 



p t p J 



M(E + M) 



It is important to note that the wave func- 
tions entering the weak current matrix element 



(3.1) are not in the rest frame in general. For 
example, in the B meson rest frame, the D** me- 
son is moving with the recoil momentum A. The 
wave function of the moving D** meson \1/ d** a 
is connected with the D** wave function in the 
rest frame ^d**o by the transformation fLU 

^.a(p)=^ c 1/2 « a X 2 « a )^«o( P ), 

(3.4) 

where R w is the Wigner rotation, Ia is the 
Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D 1 ^ 2 (R) in 
spinor representation is given by 

I °^D 1 c / q 2 (R¥j = S- 1 (p c , q )S(A)S(p), 

(3.5) 

where 



S(P) 



<P) 



2m 



1 



ap 



<P) 



is the usual Lorentz transformation matrix of the 
four-spinor. For electroweak B meson decays to 
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5-wave final mesons such a transformation con- 
tributes at first order of the I/tuq expansion, 
while for the decays to excited final mesons it 
gives a contribution already to the leading term 
due to the orthogonality of the initial and final 
meson wave functions. 

Now we can perform the heavy quark expan- 
sion for the matrix elements of B decays to ex- 
cited D mesons in the framework of our model 
and determine leading and subleading Isgur-Wise 
functions. We substitute the vertex functions 
rW and r( 2 ) given by Eqs. (O) and (O) in the 



decay matrix element (3.1) and take into account 
the wave function properties ( |3.4| ). The resulting 
structure of this matrix element is rather compli- 
cated, because it is necessary to integrate both 
over d 3 p and d 3 q. The S function in expression 



(3.2) permits us to perform one of these integra- 
tions and thus this contribution can be easily cal- 
culated. The calculation of the vertex function 
I^ 2 ) contribution is more difficult. Here, instead 
of a S function, we have a complicated struc- 
ture, containing the Qq interaction potential in 
the meson. However, we can expand this con- 
tribution in inverse powers of heavy (b, c) quark 
masses and then use the quasipotential equation 
in order to perform one of the integrations in the 
current matrix element. We carry out the heavy 
quark expansion up to first order in I/tuq. It 
is easy to see that the vertex function con- 
tributes already at the subleading order of the 
1/rriQ expansion. Then we compare the arising 
decay matrix elements with the form factor de- 



compositions ( 2.14 ) for decays to radial excita- 
tions and the corresponding ones for decays to 
orbital excitations (!(]] and determine the form 
factors. We find that, for the chosen values of our 
model parameters (the mixing coefficient of vec- 
tor and scalar confining potential e = — 1 and the 
Pauli constant k = —1), the resulting structure 
at leading and subleading order in 1/rriQ coin- 
cides with the model-independent predictions of 
HQET. This allows us to determine leading and 
subleading Isgur-Wise functions. 

4. Semileptonic decays to orbitally 
excited states 



subleading Isgur- Wise functions of semileptonic 
B decays to orbitally excited D mesons |]l6| : 
(i) B — > Diev and B — » D\ev decays 



t(w) 



1 



3 (to + l) 3 / 2 



d A p 

(2^)3 
-2e 



i>D(3/2) P 



2c„ 



Tl(w) = 
72 (w) = 



"dp ' 

A' + A 

w + 1 

w 



M D ( 3/2 ) (w - 



w - 



t(w), 

(A' + A)t(io). 



(u)B- 
, V2 



1 

D^ev and B 



1 



(4.1) 

(4.2) 
(4.3) 



D\ev decays 



3 (w + iy/ 2 
d 3 P 7 

VD(l/2) \ P 



(27T) 



A* + A 



2p 



2e„ 

^s(p), 



(w + 1) 



CM 



(4.4) 



(4.5) 



w + 1 

The contributions of all other subleading form 
factors, r)i(w) and Xi( w )i to decay matrix el- 
ements are suppressed by an additional power 
of the ratio (w — l)/(w + 1), which is equal to 
zero at w = 1 and less than 1/6 at w max = 
(1 + r 2 )/(2r). Since the main contribution to 
the decay rate comes from the values of form 
factors close to w = 1, these form factors turn 
out to be unimportant. This result is in agree- 
ment with the HQET-motivated considerations 
fic(| that the functions parametrizing the time- 
ordered products of the chromomagnetic term 
in the HQET Lagrangian with the leading order 
currents should be small. 

The arrow over d/dp in ( |4.1| ) and (4.4) indi- 
cates that the derivative acts on the wave func- 
tion of the D** meson. All the wave functions 
and meson masses have been obtained in || by 
the numerical solution of the quasipotential equa- 
tion. We use the following values for HQET pa- 
rameters A = 0.51 GeV, A' = 0.80 GeV, and 
A* = 0.89 GeV §. 

The last terms in the square brackets of the 
expressions for the leading order Isgur-Wise func- 



We get the following expressions for leading and tions t(w) (4T) and C( w ) (4.4) result from the 
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wave function transformation ( p.4| ) associated 
with the relativistic rotation of the light quark 
spin (Wigner rotation) in passing to the moving 
reference frame. These terms are numerically im- 
portant and lead to the suppression of the C form 
factor compared to r. Note that if we had applied 
a simplified non-relativistic quark model |5|, [lTj 
these important contributions would be missing. 
Neglecting further the small difference between 
the wave functions ipD(i/2) an( i V , d(3/2)> the fol- 
lowing relation between r and £ would have been 
obtained JlOfl 



ratio 



CH 



w + l 



t(w). 



(4.6) 



However, we see that this relation is violated if 
the relativistic transformation properties of the 
wave function are taken into account. At the 
point w = 1, where the initial B meson and final 
D** are at rest, we find instead the relation 



r(l)_C(l) 
V3 2 



d 3 p 



1p D *.(p)- 



P 



-ipB(p), 



(4.7) 



obtained by assuming ■0_D(3/2) — ^D(i/2) — ipD* 



The relation (4.7) coincides with the one found 
in Ref. where the Wigner rotation was also 
taken into account. 

In Table [l] we present our numerical results 
for the leading order Isgur-Wise functions r(l) 
and at zero recoil of the final D** meson, 
as well as their slopes py 2 and p\j V in com- 
parison with other model predictions [l8|, [ll]. 
|c], H, ||, [||. We see that most of the above 
approaches predict close values for the function 
r(l) and its slope p\j V while the results for ((1) 
significantly differ from one another. This dif- 
ference is a consequence of a specific treatment 
of the relativistic quark dynamics. Nonrelativis- 
tic approaches predict C(l) — {2/V3)t(1), while 
the relativistic treatment leads to (2/V3)r(l) > 
C(l). The more relativistic the light quark in the 
heavy-light meson is, the more suppressed C is 
with respect to r. 

We can now calculate the decay branching 
ratios by integrating double differential decay ra- 
tes. Our results for decay rates both in the in- 
finitely heavy quark limit and taking account of 
the first order 1/rriQ corrections as well as their 



R = 



Bt{B -► D**ev\ 



ithl/r. 



Br(£ -> D^ev)^^ 

are presented in Table |2| We see that the inclu- 
sion of 1/rriQ corrections considerably influences 
the results and for some decays their contribution 
is as important as the leading order contribution. 
This is the consequence of the vanishing of the 
leading order contribution to the decay matrix el- 
ements due to the heavy quark spin-flavour sym- 
metry at zero recoil of the final D** meson [10 , 
while nothing prevents 1 /tuq corrections to con- 
tribute to the decay matrix element at this kine- 
matical point. In fact, matrix elements at zero 
recoil are determined by the form factors fv 1 (1), 
<7+(l) and <7Vi(l), which receive non- vanishing 
contributions from first order heavy quark mass 
corrections: 



\/6/Vi(l) = -8e c (A' - A)r(l) 



(4.8) 



9+0) = -^(e c + e b )(A* -A)C(l) (4.9) 
g Vl (l) = (e c - 3efc)(A* - A)C(l). (4.10) 

Since the kinematically allowed range for these 
decays is not broad (1 < w < w max ~ 1-32), 
the contribution to the decay rate of the rather 
small 1 /rriQ corrections is substantially increased 
fll(i[| . This is confirmed by numerical calcula- 
tions. From Table ^ we see that the decay rate 
B — > D^ev, for which all contributions vanish 
at zero recoil, is only slightly increased by sub- 
leading 1/mq corrections. On the other hand, 
B — > D\ev and B — > D^ev decay rates receive 
large 1/mQ contributions. The situation is differ- 
ent for the B — > D\ev decay. Here the 1 /tuq con- 
tribution at zero recoil is not equal to zero, but it 
is suppressed by a very small factor (e c — 3sb) (see 
Eq. ( 4.10 )), which is only w 0.015 GeV -1 for our 
model parameters. As a result the B — * D\ev 
decay rate receives 1/toq contributions compa- 
rable to those for the B — > D\ev rate. The 
above discussion shows that the sharp increase of 
B — > D\ev and B — > Dq decay rates by first or- 
der 1 / mQ corrections does not signal the break- 
down of the heavy quark expansion, but is rather 
a result of the interplay of kinematical and dy- 
namical effects. Thus we have good reasons to 
expect that higher order 1/rriQ corrections will 
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Ours 


io| 


m 


m 


|21] 




[18j,^3| 


r(l) 


0.85 


0.71 


0.97 


1.14 




1.02 


0.90 


P3/2 


1.82 


1.5 


2.3 


1.9 




1.5 


1.45 


C(i) 


0.59 


0.82 


0.18 


0.82 


0.70 ±0.16 


0.44 


0.12 


Pi/2 


1.37 


1.0 


1.1 


1.4 


2.5 ± 1.0 


0.83 


0.73 



Table 1: The comparison of our results for the values of the leading Isgur-Wise functions r and £ at zero 
recoil of the final D** meson and their slopes p 2 j with other predictions. 







— * 00 


With 1/rriQ 




Experiment 


Decay 


r 


Br 


r 


Br 


R 


Br (CLEO) § Br (ALEPH) § 


B -> Diev 


1.4 


0.32 


2.7 


0.63 


1.97 


0.56 ±0.13 ±0.08 ±0.04 0.74 ±0.16 


B -> D* 2 ev 


2.1 


0.51 


2.5 


0.59 


1.16 


<0.8 <0.2 


B ->■ D\ev 


0.31 


0.073 


0.39 


0.09 


1.23 




B -> D^ev 


0.25 


0.061 


0.59 


0.14 


2.3 





Table 2: Decay rates F (in units of |V c f,/0.04| 2 x 1CT 15 GeV) and branching ratios BR (in %) for B -> D**ev 
decays in the infinitely heavy quark limit with the account of first order 1/mq corrections. R is a ratio of 
branching ratios with the account of 1/rriQ corrections to branching ratios in the infinitely heavy quark limit. 



influence these decay rates at the level of 10 - 
20%. 

In Table || we present the experimental data 
from CLEO § and ALEPH @, which are avail- 
able only for the B — > D\ev decay. For B — > 
D^ev, these experimental groups present only 
upper limits, which require the use of some ad- 
ditional assumptions about the hadronic branch- 
ing ratios of the D 2 meson. Our result for the 
branching ratio of the B — > D\ev decay with the 
inclusion of XjrriQ corrections is in good agree- 
ment with both measurements. On the other 
hand, our branching ratio for the B — > D 2 ev de- 
cay is only within the CLEO upper limit and dis- 
agrees with the ALEPH one. However, there are 
some reasons to expect that the ALEPH bound 
is too strong Q. 

Finally we test the fulfilment of the Bjorken 
sum rule IE4| in our model. This sum rule states 



P 



E 



|C<"0(i)| s 



| r M(l)p 



4 3 

m m 

(4.11) 

where p 2 is the slope of the B — » D^*'ev Isgur- 
Wise function, Q- m > and r^" 1 -' are the form fac- 
tors describing the orbitally excited states dis- 
cussed above and their radial excitations, and 
ellipses denote contributions from non-resonant 
channels. We see that the contribution of the 



lowest lying P-wave states implies the bound 
1 , IC(1)| 2 , Jr(l) 12 



P 2 > 



0.81, (4.12) 



which is in agreement with the slope p 2 = 1.02 in 
our model || and with experimental values (25). 

5. Semileptonic decays to radially ex- 
cited states 

In the case of semileptonic B decays to radially 
excited D mesons we get the following expres- 
sions for leading and subleading Isgur-Wise func- 
tions [Efl: 



£ (1) H = 



w + 1 



1/2 



(27T) 



(5.1) 

,'AW+A 1AW-A 

e«H, (5.2) 



x 1 + -- 



2 

2w-l 



3w + l 



+ - 



AW 



1 io-IAW _a 

20w±l w- 1 

1/2 



2 \w + l 
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(5.3) 

X3(W)= -80 W + 1 u,-l $ H 



AW 

4 \w + 1 



1/2 



<i 3 p 



(2tt) 3 rc 

Xfc(w) = A 



2e„ 



ID + 1 



M DW ,(w + i) 

1/2 



A)V4 0) (P), 
(5.5) 



-^ (0) /p 
^ 1)si (p)-3^ 1)sd (p) 



2e 



M D( .),(u; + 1) 



(5.6) 



where A 2 = M^ w ,(«; 2 — 1). Here we used the 
expansion for the 5*-wave meson wave function 



M 



where ip^ is the wave function in the limit tuq —>■ 
oo, ipjrf sl and V'm' 5 '' are the spin-independent 
and spin-dependent first order I/tuq corrections, 
dp = — 3 for pseudoscalar and dy = 1 for vector 
mesons. The symbol = in the expressions (5.3)- 
(5.6) for the subleading functions Xi( w ) means 



l)sd 
M 



that the corrections suppressed by an additional 
power of the ratio (w—l)/(w + l), which is equal 
to zero at w = 1 and less than 1/6 at w max , were 
neglected. Since the main contribution to the 
decay rate comes from the values of form factors 
close to w = 1, these corrections turn out to be 
unimportant. 



It is clear from the expression (5.1) that the 
leading order contribution vanishes at the point 
of zero recoil (A = 0, w — 1) of the final flW' 
meson, since the radial parts of the wave func- 
tions 'J £>(»)/ and i& b are orthogonal in the in- 
finitely heavy quark limit. The 1/toq corrections 



to the current (2.11) also do not contribute at 
this kinematical point for the same reason. The 
only nonzero contributions at w = 1 come from 
corrections to the Lagrangian 3 xi (w) , X3 i w ) an d 



Xb{w). From Eqs. ( 2.14 ) one can find for the form 



3 There are no normalization conditions for these cor- 
rections contrary to the decay to the ground state 
mesons, where the conservation of vector current requires 
their vanishing at zero recoil |12|. 



factors contributing to the decay matrix elements 
at zero recoil 

Ml) = £ C [2xi (1) + 12X3(1)] + £6X6(1), 
h Al (l) =e c [2xi(l)- 4X3(1)] +£6X6(1)- (5.7) 

Such nonvanishing contributions at zero recoil 
result from the first order 1/mg corrections to 
the wave functions (see Eq. ( |5.6| ) and the last 
terms in Eqs. ( |5.3| ), ( |5.5| )). Since the kinemati- 
cally allowed range for these decays is not broad 
( 1 < w < lUmax ~ 1.27) the relative contribution 
to the decay rate of such small 1 jvciQ corrections 
is substantially increased. Note that the terms 
£q(AW - A)£(»)(uj)/(w - 1) have the same be- 
haviour near w = 1 as the leading order contri- 
bution, in contrast to decays to the ground state 
mesons, where 1/toq corrections are sup- 
pressed with respect to the leading order contri- 
bution by the factor [w — 1) near this point (this 
result is known as Luke's theorem |l^| ). Since 
inclusion of first order heavy quark corrections 
to B decays to the ground state DM mesons re- 
sults in approximately a 10-20% increase of decay 
rates || |l3| , one could expect that the influence 
of these corrections on decay rates to radially ex- 
cited DM' mesons will be more essential. Our 
numerical analysis supports these observations. 

We can now calculate the decay branching 
ratios by integrating double differential decay 
rates. Our results for decay rates both in the 
infinitely heavy quark limit and taking account 
of the first order 1 / tuq corrections as well as their 
ratio 



Rl 



Br(g->P(*>' ei /) withl/mt 
Br(_B -> DM' e2 /) 



TTiQ — >oc 



are presented in Table |j. We find that both 
1 / toq corrections to decay rates arising from cor- 
rections to HQET Lagrangian fl5.3|)-(|5~6|), which 
do not vanish at zero recoil, and corrections to 
the current (|5.2|), fl2.ll ), vanishing at zero re- 
coil, give significant contributions. In the case of 
B — > D'ev decay both types of these corrections 
tend to increase the decay rate leading to approx- 
imately a 75% increase of the B — ► D'ev decay 
rate. On the other hand, these corrections give 
opposite contributions to the B — > D*'ev decay 
rate: the corrections to the current give a neg- 
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m Q - 


-> 00 


With 1/toq 




Decay 


r 


Br 


r 


Br 


R' 


B -> D'ev 


0.53 


0.12 


0.92 


0.22 


1.74 


B -> L>*'ei/ 


0.70 


0.17 


0.78 


0.18 


1.11 


£(£ -> D^'ev) 


1.23 


0.29 


1.70 


0.40 


1.37 



Table 3: Decay rates T (in units of |14t,/0.04| 2 x 10~ 15 GeV) and branching ratios BR (in %) for B decays 
to radially excited D^*'' mesons in the infinitely heavy quark limit and taking account of first order 1/mq 
corrections. E(_B — » D^*''ev) represent the sum over the channels. R' is a ratio of branching ratios taking 
account of 1/rriQ corrections to branching ratios in the infinitely heavy quark limit. 



ative contribution, while corrections to the La- 
grangian give a positive one of approximately the 
same value. This interplay of 1/mq corrections 
only slightly (ss 10%) increases the decay rate 
with respect to the infinitely heavy quark limit. 
As a result the branching ratio for B — ► D'ev 
decay exceeds the one for B — > D*'ev after in- 
clusion of first order 1/rriQ corrections. In the 
infinitely heavy quark mass limit we have for the 
ratio Br(B -> D'ev)/Br(B -> D*'ev) = 0.75, 
while the account of 1/iwq corrections results in 
the considerable increase of this ratio to 1.22. 

In Table || we also present the sum of the 
branching ratios over first radially excited states. 
Inclusion of 1/ tuq corrections results in approx- 
imately a 40% increase of this sum. We see that 
our model predicts that 0.40% of B meson decays 
go to the first radially excited D meson states. If 
we add this value to our prediction for B decays 
to the first orbitally excited states 1.45% fl6fl , 
we get the value of 1.85%. This result means 
that approximately 2% of B decays should go to 
higher charmed excitations. 

6. Conclusions 

In this paper we have applied the relativistic quark 
model to the consideration of semileptonic B de- 
cays to orbitally and radially excited charmed 
mesons, in the leading and subleading orders of 
the heavy quark expansion. We have found an in- 
teresting interplay of relativistic and finite heavy 
quark mass contributions. In particular, it has 
been found that the Lorentz transformation prop- 
erties of meson wave functions play an important 
role in the theoretical description of these decays. 
Thus, the Wigner rotation of the light quark spin 
gives a significant contribution already at the 



leading order of the heavy quark expansion for 
decays to orbitally excited mesons. This con- 
tribution considerably reduces the leading order 
Isgur-Wise function £ with respect to r. As a 
result, in this limit the decay rates of B — > D^ev 
and B — ► D\ev are approximately an order of 
magnitude smaller than the decay rates of B — > 
D\ev and B — > D\ev. On the other hand, inclu- 
sion of the first order 1 / tuq corrections also sub- 
stantially influences the decay rates. This large 
effect of subleading heavy quark corrections is a 
consequence of vanishing the leading order con- 
tributions to the decay matrix elements due to 
heavy quark spin-flavour symmetry at the point 
of zero recoil of the final charmed meson. How- 
ever, the subleading order contributions do not 
vanish at this kinematical point. Since the kine- 
matical range for these decays is rather small, 
the role of these corrections is considerably in- 
creased. Their account results in an approxi- 
mately twofold enhancement of the B — > D\ev, 
B — > D^ev and B — ► D'ev decay rates, while the 
B — > D^ev, B — > D\ev and B — > D*'ev rates 
are increased only slightly. The small influence of 
1 / toq corrections on the latter decay rate is the 
consequence of the additional interplay of 1/ niQ 
corrections. We thus see that these subleading 
heavy quark corrections turn out to be very im- 
portant and considerably change results in the 
infinitely heavy quark limit. For example, the ra- 
tio of branching ratios Br(£> — > D^ev) /Br(B — > 
D\ev) changes from the value of about 1.6 in 
the heavy quark limit, vhq — * oo, to the value of 
about 1 after subleading corrections are included. 
Finally we find that the semileptonic B decays to 
first orbital and radial excitations of D mesons 
amount in total to approximately 2% of the B 
decay rate. 
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